x=[1/sqrt(3) 1/sqrt(3) 1/sqrt(3)];
I=[1 0;0 1];
X=[0 1;1 0];
Y=[0 -i;i 0];
Z=[1 0;0 -1];
V=(x(1)*X+x(2)*Y+x(3)*Z);
U1=expm(-i/2*V);
U2=cos(1/2)*I-i*sin(1/2)*V;
D=(U1-I)*conj(transpose(U1-I));
D=sqrtm(D);
